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Abstract. In Part I we have developed a theory for fitting 
p-mode Fourier spectra assuming that these spectra have a 
multi-normal distribution. We showed, using Monte-Carlo 
simulations, how one can obtain p-mode parameters using 
'Maximum Likelihood Estimators'. In this article, here- 
after Part II, we show how to use the theory developed in 
Part I for fitting real data. We introduce 4 new diagnos- 
tics in helioseismology: the (ra, v) echelle diagramme, the 
cross echelle diagramme, the inter echelle diagramme, and 
the ratio cross spectrum. These diagnostics are extremely 
powerful to visualize and understand the covariance ma- 
trices of the Fourier spectra, and also to find bugs in the 
data analysis code. These diagrammes can also be used to 
derive quantitative information on the mode leakage and 
noise covariance matrices. Numerous examples using the 
LOI/SOHO and GONG data are given. 
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1. Introduction 

The physics of the solar interior is known from inversion 
of solar p-mode frequencies and splittings. These measure- 
ments are derived from fitting p-mode Fourier spectra. 
Schou (1992) was the first one to assume a multi-normal 
distribution for p-mode Fourier spectra and using a real 
leakage matrix. Following this pioneering work, Appour- 
chaux et al. (1997) (hereafter Part I), generalized the the- 
oretical background for fitting p-mode Fourier spectra to 
complex leakage matrix, and included explicitly the cor- 
relation of the noise between the Fourier spectra. Using 
Monte- Carlo simulations, we showed that our fitted pa- 
rameters were unbiased. We also studied systematic errors 
due to an imperfect knowledge of the leakage covariance 
matrix. Unfortunately, a theoretical knowledge of fitting 



Send offprint requests to: thierrya@so.estec.esa.nl 



data is not enough as only real data will teach us if our ap- 
proach is correct. Contrary to fitting p-mode power spec- 
tra, the process of fitting the Fourier spectra as described 
in Part I is rather difficult to understand and visualize. 
Schou (1992) gave a few diagnostics for understanding 
how the Fourier spectra are fitted but without showing 
an easy way to visualize the covariance matrices. 

In this paper, we show how one can easily visualize 
the mode and noise correlation matrices, and then derive 
the mode leakage matrix. In the first section, we describe 
4 new diagrammes that have various diagnostics power. 
In the second section, we describe how we use those di- 
agrammes for inferring the leakage and noise covariance 
matrices for the data of the Luminosity Oscillations Im- 
ager (LOI) on board the Solar and Heliospheric Observa- 
tory (SOHO) data (Appourchaux et al, 1997), and for the 
data of the Global Oscillations Network Group (GONG) 
(Hill et al, 1996). The LOI time series starts on 27 March 
1996 and ends on 27 March 1997 with a duty cycle greater 
than 99%. The GONG time series starts on 27 August 
1995 and ends on 22 August 1996 with a 75% duty cy- 
cle. In the last section we conclude by emphasizing the 
usefulness of these diagrammes. 



2. Diagnostics for helioseismic data analysis 

The echelle diagramme was first introduced by Grec (1981). 
It is based on the fact that the low-degree modes are es- 
sentially equidistant in frequency for a given I] the typical 
spacing for I = is 136 jaUz. The spectrum is cut into 
pieces of 136 /iHz which are stacked on top of each other. 
Since the modes are not truly equidistant in frequency, the 
echelle diagramme shows up power as distorted ridges; an 
example is given in Fig. E] for the LOI/SOHO instrument 
seeing the Sun as a star. 

Another useful diagramme was introduced by Brown 
(1985), the so-called (ra, v) diagramme which shows how 
the frequency of an I, ra mode depends upon ra. Most of- 
ten this diagramme is only shown for a single n and for 
intermediate degrees I > 10. 
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Fig. 1. Amplitude spectra echelle diagramme for 1 year 
of LOI data seeing the Sun as a star. The scale is in part- 
per- million /iHz -1 / 2 (ppm /iHz -1 / 2 ). The spacing is tuned 
for Z = 0. The I = modes are at the center, the I = 2 
modes are about 10 /iHz on the left hand side, the Z = 1 
modes are about +65 /iHz on the right hand side, the 
Z = 3 can be faintly seen at about 12 /iHz from the left 
hand side of the 1 = 1. Other modes such as Z = 4 and 
I = 5 can also be seen faintly seen at -35 /iHz and +15 
/iHz, respectively. The distortion of the ridges are due to 
sound speed gradients in the solar core. 



The purpose of these diagrammes is always to show 
an estimate of the variance of the spectra. In our case we 
also want to visualize not only the variance but also the 
covariance of the Fourier spectra. Here we briefly recall 
from Part I that the observed Fourier spectra (y) can be 
related to the individual Fourier spectra of the normal 
modes (x) by the leakage matrix C^ 1,1 ) by: 



created four new diagrammes for visualizing an estimate 
of V (M '\ all having various diagnostics power: 

1. (m, v) echelle diagramme: estimate of the diagonal el- 
ements of (Z = V) 

2. cross echelle diagramme: estimate of the off-diagonal 
elements of V (M) (Z = V) 

3. inter echelle diagramme: estimate of the off-diagonal 
elements of (1/neV) 

4. ratio cross spectrum: estimate of the ratio of the ele- 
ments of ) 

Each diagnostic is described hereafter in more detail. 



2.1. (m, v) echelle diagramme 

The (m, v) echelle diagramme is made of 21 + 1 echelle 
diagrammes of each Z,m power spectra or |^, m (^)| 2 . The 
2Z + 1 echelle diagrammes are stacked on top of each other 
to show the dependence of the mode frequency upon m. 
These diagrammes give an estimate of the diagonal of the 
covariance matrix of the observations as: 
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can be derived from the sub- matrix V^' z ) whose elements 
can be expressed as: 
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where E is the expectation, f^niy) is the profile of the 
l",m" mode, B^ 1 ) is the covariance matrix of the noise, 
and with the i//, m (^) having a mean of zero. The factor 
2 comes from the fact that the real part of V^ 1,1 ^ repre- 
sents both the covariance of the real or imaginary parts of 
the Fourier spectra (See Part I, section 3.3.2); the same 
property applied to the imaginary part of V^ 1,1 ) which rep- 
resents the covariance between the real part and the imag- 
inary part of the Fourier spectra. Equation (2j contains all 
the information that we need for visualizing an estimate 
of the real and imaginary parts of Drawing from 



where V^'^ symbolizes the estimate of V^'^. It is impor- 
tant when one makes these diagrammes to tune the spac- 
ing for the degree to study. The spacing for a given I can 
be computed from available p-mode frequencies. Since the 
spacing varies with the degree, other modes with a sig- 
nificant different spacing can be seen more like diagonal 
ridges crossing the (m, v) diagrammes; this is a powerful 
tool to identify other degrees. 

Nevertheless, the diagnostics power of the m, v echelle 
diagramme is rather limited for deriving the leakage ma- 
trix: it can be shown using Eq. (0) that the diagonal ele- 
ments of V can be expressed as: 

5 m"=l 

V&»= E \^LA 2 fU-)+B^ m (4) 

m" = — I 

As we can see with Eq. J3J, the sign information of the 
elements of C^ 1 ^ is lost; second, their magnitude being 
typically less than 0.5, the leakage elements cannot be eas- 
ily seen in the power spectra. Another kind of diagramme 
that preserves the sign of the leakage elements had to be 
devised. 



2.2. Cross echelle diagramme 

The cross echelle diagramme of an Z,m mode is made 
of 2Z + 1 echelle diagrammes of the cross spectrum of 
m and m! or yi,m(v)y* m r(i'). The 2/ + 1 real (or imagi- 
nary) parts of the cross spectra are stacked on top of each 
other to show the dependence upon m of the mode fre- 
quency. These diagrammes give an estimate of the rows 



as: 

W%% l ,(v)=y lim (v)yl m ,(y) (5) 

Of course when m = m! we get the echelle diagrammes of 
the previous section. Only I + 1 cross echelle diagrammes 
are shown as the matrix V^' z ) is hermitian by definition. 

The imaginary part of the cross spectra has some di- 
agnostic power: it represents the correlation between the 
real and imaginary parts of the Fourier spectra. When the 
leakage matrix is real, which is generally the case, there 
is no correlation between the real and imaginary parts. 
Nevertheless the imaginary part could be helpful to find 
errors in the filters applied to the images (See Part I, sec- 
tion 3.3.1). 

It can be shown that the elements of V^'^ can be ex- 
pressed as: 



V, 



(1,1) 

m,m' 



/ , u m' ,m" l -m,m" )m"\ u ! T O, 



(6) 

m" T^m' ,m 

As we can see with Eq. (0) , these diagrammes preserve the 
sign of the leakage matrix elements. In general, the cross 
spectra for m,m', representing V^'^/ carries information 

over the sign of the leakage elements C^f^, and C^'/^. 
The other additional terms expressed as product of leakage 
elements are sometimes more difficult to interpret. 

But the power of these diagrammes is not only re- 
stricted to checking the sign of the elements of the leak- 
age matrix. They are also real tools to get a first order 
estimate of the leakage matrix. We have shown in Part I, 
that there is no difference between fitting data for which 
the leakage matrix is the identity, and data for which the 
leakage matrix is not. We showed in part I, that the co- 
variance matrix of x (x = C~ 1 y) can be written, similarly 
to that of Eq. (J2J), by the sum of 2 matrices: the first one 
represents the mode covariance matrix and is diagonal, 
and the second term represents the covariance matrix of 
the noise. Therefore, applying the inverse of the leakage 
matrix to the data should, in principle, remove all artificial 
mode correlations between the Fourier spectra of x: this 
can be verified using the cross echelle diagrammes. This is 
the most powerful test for deriving the leakage matrices. 

The cross echelle diagrammes are useful to verify the 
correlation within a given degree, but other degrees are 
known to leak into the target degree, such as I — 6, 7 and 
8 into I = 1, 4 and 8, respectively. The purpose of the next 
diagramme is to assess the magnitude of these leakages. 



give an estimate of the rows (or columns) of the covariance 
matrix of the observations as: 

2V m , m /(i/) =%mMft/, m /M (?) 
Similarly as for the cross echelle diagramme, it will help to 
visualize the covariance matrix between different degrees, 
and to derive leakage elements of the full leakage matrix 
£(U ) m One can derive an equation similar to Eq. (0) for 
different degrees showing that the inter echelle diagramme 
carries information over the sign of the leakage elements 

u m,m' d,1AU u m',m ■ 

As mentioned above applying the inverse of the leakage 
matrix will help to verify to the first order that there is no 
artificial correlation due to the p modes. When different 
degrees are involved the full leakage matrix C^ 1,1 ) has to be 
used, producing diagrammes that should have no artificial 
correlation due to the p modes. 



2.4- Ratio cross spectrum 

All the previous diagrammes are helpful to understand 
and visualize the mode covariance matrices. Unfortunately, 
due to the high signal-to-noise ratio, these diagrammes 
cannot be used to evaluate the correlation of the noise in 
the Fourier spectra. In between the p modes, these corre- 
lations can be more easily visualized as we have: 
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But instead of visualizing B^ 1,1 \ we prefer to look directly 
at the correlation by computing the ratio of the cross spec- 
tra as: 
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This ratio is called the ratio cross spectrum. The ratio 
cross spectrum gives an estimate of the ratio matrix VS 1 ' 1 ) 
(See Part I) which gives a better understanding of how 
much the noise background is correlated between the Fou- 
rier spectra. Nevertheless, in the p-mode frequency range, 
the ratio cross spectrum is more difficult to interpret as 
the noise correlation is affected by the presence of the 
modes. By looking away from the modes (at high or low 
frequency) or by looking between the modes, one could 
obtain a reasonable good estimate of the noise correlation. 



2.3. Inter echelle diagramme 



3. Application to data 

3.1. (m,v) echelle diagramme for the LOI/SOHO data 

The inter echelle diagramme of an l,m mode for the de- 
gree V is made of 21' + 1 echelle diagrammes of the cross Example of these diagrammes can be seen in Figs. |2j |2|and 
spectrum of Z, m and V , m! or y^mMyf, m /(^)« The 21' + 1 0]for 1 year of LOI data for I = 1, 2 and 5, respectively, 
real part of the cross spectra are stacked on top of each It is important when one makes such diagrammes to tune 
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Fig. 2. (m, v) echelle diagramme for 1 year of LOI data for / = 1. The scale is in ppm/iHz -1 / 2 . The spacing is tuned 
for I = 1. (Left) The full diagramme centered on I = 1. The I = 3 modes are located about 15 /iHz at the left hand 
side of the 1 = 1. The I = 2 modes are on the right edge, while the I = are on the opposite side. The I = 4 modes can 
be seen around +40 /iHz. (Right) The same diagramme but enlarged around 1 = 1. The frequency shift or splitting of 
the modes due to the solar rotation can be seen: the 2 patches of power for m = — 1 and m = +1 are slightly displaced 
from each other. The artificial correlation between m = — 1 and m = +1 is not as clear. 
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Fig. 3. (ra, z/) echelle diagramme for 1 year of LOI data for I = 2. The scale is in ppm/iHz -1 / 2 . The spacing is tuned 
for I = 2. (Left) The full diagramme centered on I = 2 with a spacing of 135. /iHz. The I = modes are located about 
10 /iHz at the right hand side of the / = 2 modes. The I = 3 modes are on the right edge, while the I = 1 modes are on 
the opposite side. The I = 4 modes are easily seen at about -25 /iHz; the I = 5 starts to appear at +25 /iHz. (Right) 
The same diagramme but enlarged around I = 2 with the I = modes at the right hand side. The splitting of the 
modes is clear. Note the absence of the I = modes for m = ±1 



the spacing for the degree to study. For example, one can 
see in Fig. that the ridges of power of I = 0, 1, 2 and 3 
have different shapes than in Fig. 03 Other modes with a 
significant different spacing can be seen more like diagonal 
ridges crossing the (m, v) diagrammes; this is a powerful 
tool to identify other degrees. In Fig. El the I = 5 (m, v) 
echelle diagramme is clearly contaminated by an other 
degree, i.e. I = 8, which appears at different frequencies 
depending on rn. For I = 5, m = —5, the I = 8, m = +8 is 
quite strong; while for I = 5, m = +5, this is I = 8, m = — 8 
which shows up. This kind of 'anti'-splitting behaviour is 



typical of any aliasing degrees. It is more prominent in the 
LOI data because of the undersampling effect due to the 
large size of the LOI pixels. 



3.2. A useful detail 

Before using the other diagrammes on real data, we need 
to point out a very important property coming from the 
way the m signals are built. If the weights W^ m applied 
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Fig. 4. (m, z/) echelle diagramme for 1 year of LOI data 
for I = 5. The scale is in ppm/iHz _1//2 . The spacing is 
tuned for I = 5. Here we can clearly see the splitting 
of the I = 5 modes. Unfortunately, the I = 5 data of 
the LOI are heavily polluted by the presence of the I = 
8 modes. But since the spacing for the I = 8 modes is 
different compared with that of the / = 5 modes, the I = 8 
modes are seen as ridges going from lower left to upper 
right in each m panel; while the I = 5 modes are seen 
more like straight ridges in each m panel. It can also be 
seen that the \m\ = 5 Fourier spectra contain information 
both about m = ±5. This is the result of the large LOI 
pixel size creating oversampling. The I = 8 modes are 
also distinguished because the 'splitting' of this alias is in 
the opposite direction compared with that of the I = 5 
modes. 
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Fig. 5. Effect of a wrong sign on the cross echelle diagrammes for 1 year of LOI data for I = 2,m = —2. The scale 
is in ppm 2 /iHz _1 . The spacing is tuned for I = 2. (Left) For the real parts or Re(^ m (V))Re(^ m /(^)). (Right) For 
the imaginary parts or Im{yi^ rn {y))Im{y^ rn '{v)) without having the sign corrected. The real part of the cross echelle 
diagramme is constructed by the sum of those 2 quantities; compensation due to different signs will be misleading. 
The correlation of the real parts will always have the correct sign. This can be used for checking that the correlation 
of the imaginary parts for m = and m = 2 is different from that of the correlation of the real parts. 



Table 1. Theoretical leakage matrix of l=b for the LOI data. 
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properties as of the spherical harmonics then we have: 

W Um = Wl_ m (10) 

which means that both the Fourier spectra of +ra and 
— m can be obtained from using a single filter W/ j+m . The 
Fourier spectrum of +ra will be in the positive part of the 
frequency range, while that of — m will be in the negative 
part. This approach was first used by Rhodes et al. (1979) 
for measuring solar rotation, and mentioned theoretically 
by Appourchaux and Andersen (1990) for the case of the 
LOI. Due to the property of the Fourier transform, we 
recover, in the negative part of the frequency range, not 
the Fourier spectrum of — m but its conjugate. This fact is 
very important, if one does not take care of the sign of the 
imaginary part of — m, it will lead to very serious prob- 
lem. Needless to say that fitting data without taking into 
account this detail will have devastating effects. Figure 
shows 2 cross echelle diagrammes for the real and imagi- 
nary parts of the spectra of I = 2; the latter part had no 
sign compensation for m < 0. Obviously this important 
detail cannot be detected in the power spectra. 

3. 3. Leakage matrix measurement for a single degree 
3.3.1. LOI/SOHO data 

According to theoretical computation of the p-mode sen- 
sitivities of the LOI pixels (Appourchaux and Andersen, 
1990) and using the real shape of the LOI pixels (Appour- 
chaux and Telljohann, 1996), the leakage matrices of I = 1 
and 2 are given by: 



C(2,2) 



with: 




(11) 



= a 3 = -0.308, a 2 0.576, a 4 
—0.216. The leakage matrix of I = 5 is given in Tabled 
These leakage matrices are mean value over one year. The 
leakages vary throughout the year because the distance 
between SOHO and the Sun varies. There is no B angle 
effect as the mean B is zero over 1 year. All the leakage el- 
ements are real as the LOI filters have the same symmetry 
as the spherical harmonics. 

Figure El and |HI display the cross echelle diagramme 
for I = 1,2 and 5, respectively. From Figs. and we 
can directly verify using Eq. the sign, and sometimes 
even the magnitude of the leakage elements. The best cross 
check that our theoretical computations are correct is to 
apply the inverse of the leakage matrices to the original 
data (See Section 2.2). Figures l9l and ITUl show how the 



the LOI data; for the latter we also cleaned the original 
data from the presence of the I = modes. We shall see 
later on with the GONG data that cleaning, similar to 
that of the LOI, can be achieved not only for 2 degrees 
but also for 3 degrees (I — 1,6 and 9). 

Unfortunately, the leakage matrix of I = 5 for the LOI 
data cannot be inverted. This is the result of the pixel 
under sampling and has two dramatic consequences: the 
leakage matrix cannot be verified as for the I = 2, and 
fitting the data as described in Part I is not valid anymore 
as the leakage matrix needs to be invertible. It means that 
the 21 + 1 Fourier spectra of I = 5 are dependent. One 
way around the problem is to restrict the fitting to truly 
independent Fourier spectra. This problem is specific to 
the LOI data and applies only for I > 4. 



3.3.2. GONG data 

The leakage matrices of the GONG data have been com- 
puted by R.Howe (1996, private communication). We also 
computed similar leakage matrices using the equations de- 
veloped in Part I. The integration was made in the <fi 
plane with Omax = 65.2°, (/> ma;E = 53.5° without apodiza- 
tion. We also took into account the effect of subtracting 
the velocity of the Sun seen as a star which affects GONG 
instrument's sensitivity to modes only detected by inte- 
grated sunlight instrument (mainly the modes for which 
I + m is even). The theoretical leakage matrices of I = 1 
and 2 for GONG are also given by Eq. ([TT]) but with: 
a = -0.55, ai = -0.268, a 2 = 0.451, a 3 = -0.122, and 
«4 = —0.290. For I = 1, the leakage between m = — 1 and 
m — +1 is negative due to the subtraction of the mean 
velocity which affects these modes. For I = 2, the leak- 
age between m = — 1 and m = +1 has about the same 
value as for the LOI; these modes are not affected by the 
subtraction. 

The cross echelle diagrammes of the I — 2 GONG data 
are very close to those of the LOI (See Fig. EJ}. Figure ITT1 
shows an example of a GONG cross echelle diagramme af- 
ter having applied the inverse of the leakage matrix. It is 
clear that the p-mode correlations are removed. For I = 1 
we have used the cross echelle diagramme of m = —1 and 
m = 1 for inferring quantitatively the off-diagonal leakage 
element a. We first applied the inverse of a leakage matrix 
to the I = 1 data, and then constructed the cross echelle 
diagramme of m = — 1 and m! — 1 for these data. We then 
collapsed this diagramme by adding up all the modes with 
n = 10 — 26, and finally we corrected the collapsed dia- 
gramme from the solar noise background. The collapsed 
diagramme is shown in Fig. El (Left) for no correction 
(a = 0) and for an a of -0.53; the corrected surface of the 
collapsed diagramme as a function of a is shown on Fig. 



Fig. 6. Real part of the cross echelle diagrammes for 1 year of LOI data for I = 1. The scale is in ppm 2 /iHz -1 . The 
spacing is tuned for I = 1. (Left) For I = l,m = — 1. The lower panel represents the echelle diagramme of m = —1. 
The 2 other panels are the cross spectra with m = and m = 1. As predicted, there is no visible correlation between 
m = — 1 and m = while there is between m = — 1 and m = +1. Here we visualize the first line of the real part of the 
matrix V^ 1 ' 1 ^. (Right) For I = 1, m = 0. The bottom panel is the same as the middle panel of the left diagramme but 
with a different color scale. The middle panel is the power spectra of m = already shown in Fig. 2; the I = 6 modes 
are visible as a diagonal ridge crossing the ridge of the I = 1 modes. 
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Fig. 7. Real part of the cross echelle diagrammes for 1 
year of LOI data for I = 2. The scale is in ppm 2 /iHz -1 . 
The spacing is tuned for I = 2. (Top, left) For m = — 2. 
The middle panel (for m! = 0) has negative 'power' close 
to the frequencies of the m — — 2, indicating that as is 
negative. There is some indication that ol\ may be neg- 
ative. The upper panel (for m! = 2) has 2 patches of 
negative 'power' close to the frequencies of m = — 2 and 
m = +2, showing that is also negative. These 2 patches 
are clearly separated because of the mode splitting. (Top, 
right) For m = — 1. Here it is obvious that «2 is positive. 
(Opposite) For m = 0. There is now 2 patches of negative 
power in 2 opposite panels (m f = — 2 and m! = 2) due to 
the negativity of as. It is likely that a\ is also negative. 



Cross echelle diagram for 1=5, ™,= -5 (mean3). 
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Fig. 8. Real part of the cross echelle diagrammes for 1 year of LOI data for I = 5. The scale is in ppm 2 /iHz -1 . The 
spacing is tuned for I = 5. (Top, left) For m = — 5. Here we can see the strong leakage of m! = 5 into m = 5 and vice 
versa. The 2 patches of negative 'power' on the upper panel (or positive on the lower panel) are the m = ±5 modes 
separated by about 4-5 /iHz. (Top, right) For m = —4. Leakages between m = — 4 and m = — 2 are visible, but weaker 
than for m = —5. (Middle, left) For m = —3. Strong leakages between m = — 3 and m = +3. (Middle, right) For 
m = — 2. (Bottom, left) For m = — 1. (Bottom, right) For m = 0. The left diagrammes show more leakages compared 
with the right ones, the I + m parity makes the difference. 



Fig. 9. Real part of the cross echelle diagrammes of 1 year of LOI data for I = 1. The scale is in ppm 2 /iHz _1 . The 
inverse of the theoretical leakage matrix has been applied to the original data with a=0.474. (Left) For I = 1, m = —1. 
The artificial correlation between m = — 1 and m = +1 has been entirely removed (See Fig.E|for comparison). (Right) 
For I = 1, 77i = 0. There is no improvement as there was no correlation before. 
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, Fig. 10. Real part of the cross echelle diagrammes of 1 
year of LOI data for I = 2 after having applied the inverse 

1 of the leakage matrix. The scale is in ppm 2 /iHz _1 . The 
spacing is tuned for I = 2. (Top, left) For m = —2. (Top, 

2 right) For m = — 1. (Opposite) For m = 0. Using the 
inverse of the leakage matrix we have removed all artificial 

3 correlations (See Fig.[3for comparison). The spectra have 
been cleaned both from the m leaks and from I = modes 
because we used the leakage matrix of I = 2 and 1 = 0. 



H2I (Right). When the corrected surface is close to 0, there 
is no artificial correlation remaining. 

Cleaning the data from artificial correlations has also 
been done in a different way by Toutain et al. (1998). Us- 



ing Singular Value Decomposition, they recomputed pixel 
filters for the MDI/LOI proxy so as to remove the m leaks 
and the other aliasing degrees. Here we have shown, that 
data cleaning is also possible without having the pixel time 



technique is more useful as one has to reduce a smaller 
amount of data. 



3.4- Leakage matrix measurement for many degrees 
3.4.1. LOI/SOHO data 

Figures IT3l and PHI shows the inter echelle diagramme for 
I = 1, V = 6, and I = 5, V = 8, respectively. These di- 
agrammes are more difficult to interpret, but using spe- 
cific spacing for I or l\ one can unambiguously identify 
which degree contributes to the correlation. For the LOI 
data, the I — 5 and I = 8 (similarly I = 4 and I = 7) 
are contaminated by each other. This is again due to the 
spatial undersampling which prevents the LOI data to be 
cleaned from aliasing degrees making the leakage matrix 
not invertible. Fortunately, the I = 1 LOI data are far 
less affected by the presence of I = 6 compared with that 
of GONG (See in the next section); this is the result of 
an effective apodization function (limb darkening) which 
creates a narrower spatial response than that of GONG 
(line-of-sight projection). 

Figure El nas already shown that we could clean the 
I = 2 LOI data from the I = mode. This was done 
by taking into account the leakage matrix C^ ' 2 ). Unfortu- 
nately as mentioned for the I = 5, the LOI data cannot be 
cleaned from aliasing degrees because the leakage matrice 
C (i,6\ and C^ 5 ' 8 ) are singular. 



3.4.2. GONG data 

Figures El an d El show the inter echelle diagramme for 
I = 1,/' = 6, and for I = 1,Z' = 9. The correlations are 
clearly visible. If they are not taken into account they are 
likely to bias the frequency and splittings of the I = 1 
(See Rabello-Soares and Appourchaux, 1998, in prepara- 
tion). By using the inverse full leakage matrix of I = 1,6 
and 9 (i.e. C^ 1 ' 6 ' 9 ^ , one can clean the data from these spuri- 
ous correlations arising both from the aliasing degrees and 
from the m leaks. Figures IT6l and IT8l shows the 'cleaned' 
inter echelle diagramme which should be compared with 
Figs. El and El respectively. The correlations have been 
almost entirely removed. Some correlations between the 
I = l,m = — 1 modes and the I = 6, 9 modes are still 
present. This is due to the fact that the computation of 
the leakage for I = l,m = ±1 are very sensitive to the 
subtraction of the velocity of the Sun seen as a star. This 
is not the case for the I = 1, m = modes as they are in- 
sensitive to this effect. Nevertheless, our imperfect knowl- 
edge of the leakage matrix can be adjusted in order to 
remove the residual correlations. It helped to reduce the 



cies, mainly where the I = 6 start to cross over the I = 1 
(Rabello-Soares and Appourchaux 1998, in preparation). 

It is possible to apply this cleaning technique to higher 
degree modes (I « 50 — 100. In principle, this is feasible, 
although the plethora of data involved may be fairly sub- 
stantial. For higher degrees, the ridges of the modes in 
the (m, v) echelle diagrammes are almost parallel to each 
other. The idea would be not to use the full Fourier spec- 
tra, as we did for the low degree GONG data, but to use 
only a small frequency range of about ± 30 /iHz around 
the target degree. By doing so, we will not only clean the 
data from the aliasing degrees but also from the m leaks. 
In this case, the mode covariance matrix is diagonal. If we 
assume, wrongly, that the noise covariance matrix is diag- 
onal, the statistics of each cleaned spectrum, for high de- 
gree modes, could be approximated by a \ 2 with 2 degree 
of freedom. Neglecting the off-diagonal elements of the 
noise covariance matrix will lead to larger mode linewidth, 
thereby producing underestimated splittings; this system- 
atic bias will decrease as the degree. This approximation 
will produce much less systematic errors than the approx- 
imation used by Hill et al. (1996) for the original data. 
The use of the suggested approximation for the GONG 
data may help, at the same time, to reduce systematic 
errors, and to increase computing speed for higher degree 
modes. The degree where one needs to switch between 
this approximation and the correct analysis needs to be 
determined. 



3.5. Noise covariance matrix measurement 
3.5.1. LOI/SOHO data 

The leakage and ratio matrices have similar properties 
(See Part I). For example, for a given I there is no cor- 
relation between the m for which m + m' is odd for ei- 
ther matrix. This lack of correlation can be seen in Fig. 
El f° r the ratio cross spectra of the LOI data. The ratio 
matrix is also useful if one wants to reduce the number of 
noise parameters to be fitted. This is useful when the noise 
background, mainly of solar origin, varies slowly over the 
p-mode range. In this case the ratio can be assumed to be 
constant over the p-mode range. For example for I = 1 we 
can fit 2 noise parameters instead of 3. When the noise 
varies in the p-mode range, it is more advisable to fit as 
many noise parameters as required. 

Figure [201 shows the ratio cross spectra of I = 2 for the 
LOI data. As this is commonly the case for the LOI, the 
ratios are independent of frequency. In addition, as out- 
lined in Part I, the ratio matrix is very close to the leakage 
matrix. Using the values of a given for the LOI in Eq. (|11[) 
and Fig. [201 one can see that this is the case for the LOI. 
The ratio matrix has even the same symmetry property as 
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Fig. 11. Real part of the cross echelle diagrammes for 
360 days of GONG data for I = 2 after having applied the 
inverse of the leakage matrix. The scale is in m 2 s -2 Hz -1 . 
The spacing is tuned for I = 2. (Top, left) For m = — 2. 
(Top, right) For m = — 1. (Opposite) For m = 0. As for 
the LOI data the artificial correlations are removed. 
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Fig. 12. (Left) Collapsed diagramme of the real part of the cross spectrum of m = —1 and m = +1 for the GONG 
data after having applied the inverse of the leakage matrix, (top) no correction, (bottom) a=-0.53. (Right) Surface of 
the collapsed diagramme as a function of a, the surface is for a=-0.53. 



the leakage matrix. In the case of the LOI, we sometimes 
use the independence of the ratio with frequency to reduce 
the number of free parameters. It is less straightforward 
to measure the noise correlation for the GONG data than 
for the LOI data. We recommend to measure it in between 
the p modes because that is what the fitting routines will 
determine. 



3.5.2. GONG data 

Figures l2ll and l22l show the ratio cross spectra of I = 2 for 
the GONG data. The ratio matrices (as the leakage ma- 
trices) are not symmetrical. Apart from the 1 = 1, these 
matrices tend to be very close to those of the leakage ma- 
trices (see previous section). It is also clear that the ra- 
tios depend upon frequency, probably due to the effect 
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Fig. 13. Real part of the inter echelle diagrammes for 1 year of LOI data for Z = 1 and V = 6. The scale is in 
ppm 2 /iHz _1 . The spacing is tuned for V = 6. (Left) For Z = l,m = —1. Given the structure of the ridges, it means 
that the Z = 1 modes leak into the even m modes of V = 6, while the V = 6 modes do not leak into the 1 = 1. (Right) 
For Z = 1, m = 0. Again using the structure of the ridges, the odd m! of V = 6 leaks weakly into Z = 1, m = 0. 



of mesogranulation affecting the spatial properties of the 
noise with frequency. In this case, we cannot reduce the 
number of parameters to be fitted (see Rabello-Soares and 
Appourchaux 1998, in preparation) as we do for the LOI. 



4. Summary and Conclusion 

For helping the understanding of fitting p-mode Fourier 
spectra, we have devised 4 new helioseismic diagnostics: 

1. the (ra, v) echelle diagramme helps visualizing: 

— which degree can be detected, 

— the mode splitting and 

— spurious degrees. 

2. the cross echelle diagramme helps verifying: 

— the sign of the imaginary parts of the Fourier spec- 
tra for m < 0, 

— the sign of the elements of the leakage matrix of 
the degree Z, 

— to the first order the theoretical leakage matrix of 
the degree Z by applying the inverse of this matrix 
to the data (i.e. cleaning the data from the m leaks) 
and 

— the cleaning of the data from m leaks. 

3. the inter echelle diagramme helps verifying: 

— the sign of the elements of the full leakage matrix 
for degrees Z, Z', 

— to the first order the theoretical leakage matrix of 
the degree Z, V by applying the inverse this matrix 
to the data and 

— the cleaning of the data from aliasing degrees. 

4. the ratio cross spectrum helps deriving : 

— the amplitude of the noise correlation between the 
Fourier spectra and 



— the frequency dependence of the correlations. 

These steps will help to evaluate the matrices of the leak- 
age, mode covariance, ratio and noise covariance directly 
from observations. The fitting of the p modes will be con- 
siderably eased by verifying that the theoretical knowledge 
of these matrices is correct. 

These steps have been used both on data for which 
we design the spatial filtering (LOI instrument), and on 
data for which we did not design the filtering (GONG). As 
a matter of fact these diagnostics are so powerful that a 
theoretical knowledge of the various matrices is not always 
necessary to understand the data. Very often these diag- 
nostics can also be used to find bugs in the data analysis 
routines (Fig. EJ). Nevertheless, it is advisable to know to 
the first order the leakage and ratio matrices for speeding 
up the analysis process. 

Last but not least, the knowledge of the leakage ma- 
trices can be used for cleaning the data from m leaks 
and from undesired aliasing degrees. This cleaning can be 
performed either when the pixel time series are available 
(MDI data, Toutain et al. , 1998) or more simply when 
only the Fourier spectra are available (LOI, GONG). The 
cleaning has very useful application for the GONG veloc- 
ity data for removing aliasing degrees from Z = 1, 6 and 
9 for instance, and for inferring better Z = 1 splittings 
(Rabello-Soares and Appourchaux 1998, in preparation). 
The cleaning is somewhat easier with Fourier spectra as 
one does not need to reduce large amount of image time 
series. It is advisable that, in the near future, this cleaning 
technique be used for higher degree modes. This will hope- 
fully provide helioseismology with frequency and splitting 
data having much less systematic errors than before. 
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Fig. 14. Real part of the inter echelle diagrammes for 1 year of LOI data for / = 5 and V = 8. The scale is in 
ppm 2 /iHz -1 . The spacing is tuned for V = 5. (Top, left) For I = 5, m = —5. (Top, right) For I = 5, m = —4. (Middle, 
left) For / = 5,m = -3. (Middle, right) For / = 5, m = -2. (Bottom, left) For I = 5,m = -1. (Bottom, right) For 
I = 5, m = 0. The Z = 5 ridges are more like vertical lines while those of the I = 8 are more like diagonal lines going 
across from left to right. 
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Fig. 15. Real part of the inter echelle diagrammes of 1 year of GONG data for I = 1 and V = 6. The scale is in 
m 2 s _2 Hz _1 . The spacing is tuned for V = 6. (Left) For I = l,m = —1. The strongest correlations are for m! = ±6 
due to V = 6. (Right) For Z = 1, m = 0. The I = 6 modes are strong in all the panels. For m! = ±1 other modes can 
be seen as an 'asterisk' ridge: the ridge going from bottom right to upper left represents the I = 1 modes; the ridge 
going from bottom left to upper right represents the I = 9 modes; the vertical ridge represents the I = 6 modes. At 
the far left, the ridge parallel to that of the I = 1 represents the I = 3 modes. 
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Fig. 16. Real part of the inter echelle diagrammes of 1 year of GONG data for I = 1 and V = 6 after having applied 
the inverse of the leakage matrix C^ ly6,9 \ The scale is in m 2 s _2 Hz _1 . The spacing is tuned for V = 6. (Left) For 
I = 1, m = —1. The faint yellow oblique ridges are due to the I = 1 modes. (Right) For I = 1, m = 0. The correlations 
due to the I = 1 and 6 modes are almost entirely removed. 



Fig. 17. Real part of the inter echelle diagrammes of 1 year of GONG data for I = 1 and V = 9. The scale is in 
cm 2 s -2 /iHz -1 . The spacing is tuned for I — 9. (Left) For I = 1, m = — 1. The strongest correlations are for m! = ±9 
due to the V = 9 modes. (Right) For I = l,m = 0. The strongest correlations are for ra' = 0, ±1 mainly due to the 
V = 9 modes. A faint oblique ridge can be seen for m! = due to the I = 1 modes. 



Inter Echelle diagram for 1=1, m,= -l (gongtri), and 1=9 Inter Echelle diagram for 1=1, ^=0 (gongtri), and 1=9 
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Fig. 18. Real part of the inter echelle diagrammes of 1 year of GONG data for / = 1 and V = 9 after having applied 
the inverse of the leakage matrix C^ 1,6y9 \ The scale is in cm 2 s -2 /iHz -1 . The spacing is tuned for V = 9. (Left) For 
I = l,m = —1. The correlations due I = 9 are not correctly compensated for. (Right) For I = l,m = 0. The I = 3 
correlation are also enhanced; they can be seen as oblique ridges at the left hand side of the diagrammes. 
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Fig. 19. Real part of the ratio cross spectra for Z = 1. (Left) For m = — 1 and m'=0. There is very little correlation, 
less than -0.05. The modes also are not seen. (Right) For m = — 1 and m = 1. The noise correlation is about 0.55. The 
modes are visible with various positive and negative correlation. 
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Fig. 20. Real part of the ratio cross spectra of / = 2 for the LOI data. (Top, left) For m = — 2 and m! = 0. The 
noise correlation is about -0.45. (Top, right) For m = — 2 and m = 2. The noise correlation is about -0.20. (Bottom, 
left) For m — — 1 and m! = 1. The noise correlation is about 0.60. (Bottom, right) For m = and m! — 2. The noise 
correlation is about -0.35 same as for m = —2 and m! = 0. 
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Fig. 21. Real part of the ratio cross spectra of Z = 1 for the GONG data. (Left) For m = — 1 and m'=0. As for the 
LOI, there is very little correlation, less than -0.05. The modes also are not seen. (Right) For m = — 1 and m = 1. The 
noise correlation is -0.77. The correlation is negative because of the way the GONG data are reduced: the full-disk 
integrated velocity has been removed from the signal. 
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Fig. 22. Real part of the ratio cross spectra of / = 2 for the GONG data. (Top, left) For m = — 2 and m! = 0. The 
noise correlation is rather small and about -0.1. (Top, right) For m = — 2 and m = 2. The noise correlation is frequency 
dependent, with a typical value at 3000 /iHz of about -0.3. (Bottom, left) For m — — 1 and m' = 1. The typical value 
is 0.4. (Bottom, right) For m = and m! = 2. The typical value is -0.3 



